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A.  INTRODUCTION 


Leaky  wave  antennas  for  the  millimeter  wave  range  face  two  main 
problems.  The  first  relates  to  the  small  wavelengths  involved,  which 
require  small  waveguide  dimensions  and  pose  fabrication  difficulties. 
The  second  problem  is  higher  metal  loss;  for  antennas  which  are  many 
wavelengths  long,  the  leakage  (which  results  in  radiation)  may  compete 
with  the  intrinsic  waveguide  loss,  and  the  antenna  design  may  be 
adversely  affected.  In  this  report,  we  present  an  analysis  for  a 
new  type  of  leaky  wave  antenna  for  millimeter  waves  that  addresses  both 
of  these  problems.  This  new  antenna  overcomes  the  first  of  these  prob¬ 
lems  by  employing  a  longitudinally  continuous  aperture,  and  the  second 
by  basing  the  antenna  on  a  low-loss  waveguide,  the  groove  guide. 

The  groove  guide  was  one  of  several  waveguiding  structures  proposed 
for  millimeter  wave  use  about  20  years  ago  in  order  to  overcome  the 
higher  attenuation  occurring  at  these  higher  frequencies.  Although 
these  new  low-loss  waveguides  were  introduced  many  years  ago,  they  were 
not  pursued  beyond  some  initial  basic  studies  because  they  were  not  yet 
needed.  Now,  millimeter  waves  are  again  becoming  important,  and  atten¬ 
tion  is  again  being  paid  to  new  types  of  waveguides. 

The  open  groove  guide  is  shown  in  Fig.  1,  and  an  indication  of  the 
dominant  mode  electric  field  lines  present  in  its  cross  section  is  given 
in  Fig.  2(a).  One  should  first  note  that  the  structure  resembles  that 
of  rectangular  waveguide  with  most  of  its  top  and  bottom  walls  removed. 
The  groove  guide  can  therefore  be  excited  by  providing  a  smooth  tapered 
transition  between  it  and  a  feed  rectangular  waveguide.  Furthermore,  if 
symmetry  is  maintained,  many  components  can  be  designed  for  groove  guide 
which  are  analogues  of  those  in  rectangular  guide. 

With  respect  to  the  low-loss  nature  of  groove  guide,  we  should 
recall  that  when  the  electric  field  is  parallel  to  the  metal  walls  the 
attenuation  associated  with  those  walls  decreases  as  the  frequency  is 
increased;  conversely,  the  attenuation  increases  with  increasing  fre¬ 
quency  when  the  electric  field  is  perpendicular  to  the  walls.  Since  In 
groove  guide  the  electric  field  is  seen  to  be  mostly  parallel  to  the 
walls,  its  overall  attenuation  at  higher  frequencies  is  much  lower  than 
that  of  rectangular  waveguide,  where  most  of  the  field  is  perpendicular 
to  the  top  and  bottom  walls. 


Fig.  1  The  open  groove  guide,  comprised  of  two  parallel  metal  plates  whose 
central  regions  are  grooved  outwards. 


(a)  (b) 

Fig.  2  The  electric  field  of  the  dominant  mode  in  open  groove  guide,  (a)  A 

sketch  of  the  electric  field  lines  in  the  cross  section,  (b)  An  approximate 
plot  of  the  vertical  component  Ey  as  a  function  of  vertical  position  y, 
showing  that  the  guided  mode  is  hound  transversely  to  the  central  grooved 
region. 
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The  greater  width  in  the  middle,  or  central,  region  was  shown  by  T. 
Nakahara  [1,2],  the  inventor  of  groove  guide,  to  serve  as  the  mechanism 
that  confines  the  field  in  the  vertical  direction,  much  as  the  dielec¬ 
tric  central  region  does  in  H  guide.  The  field  thus  decays  exponen¬ 
tially  away  from  the  central  region  in  the  narrower  regions  above  and 
below,  as  shown  in  Fig.  2(b).  If  the  narrower  regions  are  sufficiently 
long,  it  does  not  matter  if  they  remain  open  or  are  closed  off  at  the 
ends.  - 


Work  on  the  groove  guide  progressed  in  Japan  [2,3]  and  in  the 
United  States  [4-6]  until  the  middle  1960's,  but  then  stopped  and  was 
later  revived  and  developed  further  by  D.  J.  Harris  and  his  colleagues 
[7,8]  in  Wales.  Their  recent  work  is  mainly  experimental,  being  associ¬ 
ated  with  components  for  groove  guide.  Other  recent  work  includes  that 
of  J.  Meissner,  who  investigated  the  coupling  [9]  between  groove  guides, 
and  the  radiation  losses  of  bends  [10].  To  our  knowledge,  the  present 
study  represents  the  first  contribution  to  antennas  based  on  the  groove 
guide. 

A  presentation  of  the  new  antenna  structure  and  a  brief  summary  of 
its  properties  were  made  recently  in  a  pair  of  letters  [11,12].  In  this 
comprehensive  report,  we  present  a  complete  theoretical  analysis 
together  with  numerical  results  and  systematic  design  considerations. 

The  actual  antenna  structure  and  its  operating  principle  are 
described  in  Section  B.  The  basic  principle  is  that  of  introducing  asym¬ 
metry  in  the  cross  section  so  as  to  convert  an  initially  bound  mode  into 
a  leaky  mode.  Although  it  is  applied  here  to  a  specific  structure,  this 
basic  principle  is  quite  general,  and  other  asymmetry  mechanisms  can  be 
employed  to  perform  the  same  function,  both  in  groove  guide  and  in  other 
open  waveguides  in  which  the  basic  mode  is  bound. 

The  properties  of  the  leaky  mode  (which  provides  the  radiating 
aperture  for  the  antenna)  are  obtained  from  an  accurate  transverse 
equivalent  network  the  constituents  of  which  are  derived  in  Section  C. 
In  these  derivations  it  is  necessary  first  to  identify  the  correct 
transverse  modes,  and  then  to  deduce  expressions  for  the  network  consti¬ 
tuents  which  are  accurate  but  also  simple  in  form,  so  that  the  resulting 
dispersion  relation  permits  numerical  calculations  without  undue  com¬ 
plexity.  We  were  fortunate  in  being  able  to  derive  all  the  network 
parameters  in  closed  form.  The  most  difficult  constituent  was  that  for 
the  added  asymmetric  longitudinal  strip  that  couples  the  leaking 
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transverse  mode  to  the  bound  transverse  mode.  Its  equivalent  network  is 
derived  by  employing  small  obstacle  theory  in  a  multimode  context,  where 
one  of  the  modes  is  below  cutoff,  thereby  requiring  a  modification  in 
the  expressions  available  in  the  literature. 

From  the  complete  transverse  equivalent  network,  a  dispersion  rela¬ 
tion  is  obtained  which  is  in  closed  form.  The  dispersion  equation  and 
the  numerical  procedure  employed  to  obtain  specific  numerical  results 
are  discussed  in  Section  D. 

In  Section  E  we  present  design  considerations  for  the  leaky  wave 
antenna  together  with  a  variety  of  numerical  results.  These  design 
aspects  include  optimization  considerations  for  the  leakage  constant, 
and  the  relations  between  the  features  of  the  radiation  pattern  and  the 
phase  and  leakage  constants.  The  antenna  is  seen  to  be  capable  of  sys¬ 
tematic  design  and  to  be  versatile  with  respect  to  the  scan  angle  and 
the  beam  width  of  the  radiation. 


B.  THE  OPERATING  PRINCIPLE  OF  THE  NEW  ANTENNA 

The  new  leaky  wave  antenna  is  shown  in  Fig.  3.  The  basic  differ¬ 
ence  between  the  structures  in  Figs.  1  and  3  is  that  in  Fig.  3  a  con¬ 
tinuous  metal  strip  of  narrow  width  has  been  added  to  the  guide  in  asym¬ 
metrical  .  fashion.  Without  this  strip,  the  field  of  the  basic  mode  of 
the  symmetrical  groove  waveguide  is  evanescent  vertically,  so  that  the 
field  has  decayed  to  negligible  values  as  it  reaches  the  open  upper  end. 
The  function  of  the  asymmetrically  placed  metal  strip  is  to  produce  some 
amount  of  net  horizontal  electric  field,  which  in  turn  sets  up  a  mode 
akin  to  a  TEM  mode  between  parallel  plates.  The  field  of  that  mode  pro¬ 
pagates  all  the  way  to  the  top  of  the  waveguide,  where  it  leaks  away. 
It  is  now  necessary  to  close  up  the  bottom  of  the  waveguide,  as  seen  in 
Fig.  3,  to  prevent  radiation  from  the  bottom,  and  (nonelectrically)  to 
hold  the  structure  together.  Of  course,  the  upper  walls  could  end  as 

shown  in  Fig.  3  or  they  could  attach  to  a  ground  plane. 


3  Cross  section  of  the  new  leaky  wave  antenna,  where  leakage  is  produced 
by  the  introduction  of  an  asymmetric  continuous  metal  strip  of  width  6. 


The  structure  in  Fig.  3  now  yields  a  leaky  wave  line-source  antenna 
of  simple  configuration,  since  it  is  longitudinally  continuous.  The 
value  of  the  phase  constant  6  of  the  leaky  wave  is  governed  primarily  by 
the  properties  of  the  original  unperturbed  groove  guide,  and  the  value 
of  the  attenuation,  or  leakage,  constant  a  is  determined  by  the  width 
and  location  of  the  perturbing  strip.  As  discussed  in  Section  E,  stand¬ 
ing  waves  associated  with  the  added  mode  created  by  the  asymmetry  intro¬ 
duce  additional  problems  in  the  design,  but  these  problems  can  be  over¬ 
come  by  certain  optimization  considerations. 

As  mentioned  in  the  Introduction,  the  asymmetry  can  be  produced  in 
other  ways,  resulting  in  different  actual  structures  but  with  antennas 
that  perform  similarly.  This  particular  structure  was  chosen  because  it 
is  amenable  to  accurate  analysis  in  a  simple  way,  as  explained  in  detail 
in  the  next  section. 


C.  THE  TRANSVERSE  EQUIVALENT  NETWORK 


» 


The  complete  transverse  equivalent  network  for  the  antenna  struc¬ 
ture  depicted  in  Fig.  3  is  derived  below  by  viewing  it  as  an  open 
waveguide  with  a  complex  propagation  wavenumber  8  -  jo  in  the  z  direc¬ 
tion.  The  dispersion  relation  for  8  and  a  is  obtained  in  Section  D  by 
applying  the  transverse  resonance  condition  to  this  transverse 
equivalent  network. 

In  this  transverse  equivalent  network,  the  uniform  regions  of  the 
open  waveguide's  cross  section  are  represented  by  transmission  lines, 
and  junctions  or  other  discontinuities  are  represented  by  lumped  ele¬ 
ments.  The  first  step  in  the  derivation  of  the  network  is  therefore  the 
identification  of  the  correct  transverse  modes  which  the  transmission 
lines  represent.  The  lumped  elements  must  then  be  recognized  or  derived 
corresponding  to  the  various  discontinuities  present  in  the  cross  sec¬ 
tion. 


As  discussed  in  Section  B,  the  antenna  is  created  by  introducing 
into  the  basic  nonradiating  groove  guide  an  additional  longitudinally 
continuous  metal  strip.  The  strip  is  introduced  in  asymmetrical  fashion 
so  as  to  produce  some  net  horizontal  electric  field.  In  effect,  a 
new  transverse  mode  is  created  thereby,  which  propagates  (at  an  angle) 
all  the  way  to  the  top  of  the  waveguide,  where  it  leaks  away,  thus 
transforming  the  initially  bound  longitudinal  mode  into  a  leaky  one. 
The  transverse  equivalent  network  must  therefore  be  based  on  these  two 
transverse  modes,  which  propagate  in  the  y  direction  and  are  coupled 
together  by  the  added  strip.  The  coupled  transverse  modes  then  combine, 
of  course,  to  produce  a  net  TE  longitudinal  mode  (in  the  z  direction) 
with  a  complex  propagation  constant  8  -  ja. 

The  basic  form  of  the  transverse  equivalent  network  is  shown  in 
Fig.  M,  where  the  network  has  been  placed  on  its  side  for  clarity.  In 
the  network,  the  i-1  transmission  lines  represent  the  original  mode  (the 
dominant  nonradiating  mode)  with  a  half  sine  variation  in  the  x  direc¬ 
tion  in  Fig.  3i  and  the  i=0  transmission  lines  represent  the  new  mode 
which  has  no  variation  with  x.  The  characteristic  admittances  Y  and  the 
transverse  wavenumbers  ky  of  these  transmission  lines  must  be  properly 
identified  and  related,  and  expressions  for  the  lumped  elements 
representing  the  various  discontinuities  in  the  cross  section  must  be 
derived.  The  proper  mode  identifications  and  the  mode  functions 
corresponding  to  the  transmission  lines  are  discussed  in  the  next 


i 
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subsection.  In  subsequent  subsections,  expressions  are  presented  or 
derived  for  the  step  junctions,  the  asymmetrical  coupling  strip,  and  the 
radiating  open  end.  In  the  final  subsection  of  Section  C,  the  various 
constituents  are  combined  to  produce  the  complete  transverse  equivalent 
network. 

1 .  The  Transverse  Modes 

To  properly  characterize  the  transmission  lines  in  the  transverse 
equivalent  network  in  Fig.  H,  we  must  Identify  the  correct  modes  in  the 
y  direction.  We  first  note  that  with  respect  to  the  z  (longitudinal) 
direction  the  overall  guided  mode  is  a  TE  (or  H)  mode;  that  is,  there 
exists  only  a  component  of  H  in  the  z  direction.  This  result  is  to  be 
expected  since  the  groove  guide  consists  of  a  perfectly  conducting  outer 
structure  filled  with  only  a  single  dielectric  material  (air).  In  the  y 
direction,  however,  there  exist  both  Ey  and  Hy  components,  so  that  the 
mode  is  hybrid  in  that  direction. 

Since  the  groove  guide  is  uniform  in  the  z  direction,  and  its  field 

has  only  an  H  component,  a  hybrid  mode  in  the  y  direction  is  seen  to  be 

what  is  called  by  some  an  H-type  mode  with  respect  to  the  z  direction, 

and  by  others  an  LSE  mode  with  respect  to  the  z  direction.  We  prefer 

the  former  notation,  and  we  shall  designate  a  hybrid  mode  in  the  y 
( z ) 

direction  as  an  H  -type  mode  .  Altschuler  and  Goldstone  [13]  discuss 
such  modes  in  detail  and  present  the  field  components  for  them  and  the 
characteristic  admittances  for  transmission  lines  representative  of 
them. 

For  the  transmission  lines  in  Fig.  H  we  therefore  require  the  mode 

(z) 

functions  and  transmission  line  properties  of  an  H  -type  mode  in 
parallel-plate  guide,  which  propagates  in  the  y  direction  and  is  hybrid 
in  that  direction,  but  has  only  an  Hz  component  in  the  z  direction.  The 
coordinate  system  we  employ  is  that  given  in  Fig.  3,  but  it  differs  by  a 
rotation  from  the  one  employed  in  reference  [13].  In  our  transverse 
equivalent  network,  we  are  concerned  with  two  different  -type  modes 

in  the  y  direction,  the  i  -  1  and  the  i  -  0  modes.  We  recognize  that  in 
the  sum-of-squares  relation 

2  2  2  2 

k  -kx+ky+kz  (1) 

the  longitudinal  propagation  wavenumber  k  (-0-jo)  is  the  same  for  both 

z 

modes,  and  the  unknown  k  values  are  designated  as  k  ,  and  k  _.  The  k 

y  yi  yo  x 
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values  are  k  .  -  —  and  k  «  0,  consistent  with  their  known  variations 

XI  cl  xo 

with  x.  The  characteristic  admittance  for  either  of  the  modes  is  given 
by  [13]  as 

,2,2 

k  "kz  1 


where  ky  is  of  course  the  propagation  wavenumber  of  the  transverse 
transmission  line. 

The  components  of  the  electric  and  magnetic  fields  (of  these 
transverse  modes)  that  are  transverse  to  the  y  direction  are  given  by 
the  usual  expressions 

E^x.y.z)  -  V(y)£(x,z)  (3) 

Ht(x,y,z)  -  I(y)h(x,z)  (4) 

where  the  voltage  V(y)  and  current  I(y)  satisfy  the  transmission  line 
equations  in  the  y  direction,  and  £  and  h  are  mode  functions  dependent 
only  on  the  cross  section  to  y.  Consistent  with  the  specification  (2) 
for  characteristic  admittance,  the  electric  and  magnetic  field  vector 
mode  functions,  £  and  £,  satisfy  the  orthonormality  condition 

/  £  «  y^ •£*  dS  -  1 
s 

where  the  integration  is  performed  over  the  cross  section  normal  to  ,y, 
and  we  have 

e  (x,z)  -  -h  (x,z) 

X  z 

The  z  dependence  everywhere  is  understood  to  be  exp(-jkzz). 

For  the  parallel-plate  region,  using  (5),  we  find  for  the  i-1  mode 
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We  will  need  the  hyl  component  of  the  mode  function  in  the  derivation  of 
the  expression  for  the  added  coupling  strip.  That  component  is  obtained 
from 


hy(x,z) 


yielding 


-Jky  3hz(x,z) 
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hy1(x,z) 
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Correspondingly, 

ex0(x.z)  - 

Vx,2)  ‘  |I 


for  the  i-0  mode,  we  have 


-jkzz 

e  -  -hZ0(x'z) 
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or 


hyO(X,Z) 


(10) 
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from  k  -  kyQ  +  kz  ,  since  k^  -  0.  The  fields  simplify  considerably 
for  the  i«0  mode.  As  seen,  the  mode  functions  do  not  depend  on  x,  and 
the  mode  is  in  fact  not  hybrid  but  a  TE  mode  in  the  y  direction,  having 
E  ,H  and  H  components  only.  As  a  result,  expression  (2)  for  the 

«  /  Z 

characteristic  admittance  reduces  to 
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(11) 
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The  following  complete  summary  of  the  field  components  for  the  i-1 
and  i-0  modes  is  obtained  on  use  of  (3)  and  (4)  together  with  the  mode 
functions  above  and  the  remaining  ones  obtained  in  the  same  fashion;  g 
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0 


E  (x.y.z) 
Hx(x,y,z) 


Ji,(y) 


,2,2 

k  -k 

z 


pi 

\]  a  a 


,  ‘^z2 

cos  —  x  e 
a 


Hy(x,y,z) 


H  (x.y.z) 
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We  should  note  that  the  field  ratios 


are  equal  to  the  characteristic  impedance  of  the  longitudinal  TE  mode 
(in  the  z  direction),  as  expected.  The  specific  values  of  the  V  and  I 
terms  depend  on  the  location  in  the  cros3  section,  and  are  determined  by 
the  dispersion  relation. 


2.  The  Step  Junctions 

Two  identical  step  junctions  appear  in  the  cross  section  of  the 
antenna  shown  in  Fig.  3,  but  the  effect  of  a  step  junction  on  the  i-1 
and  the  i-0  transverse  modes  is  different.  We  should  also  recognize 
that  these  step  junctions  do  not  couple  the  i-1  and  i-0  modes.  We  thus 
require  separate  expressions  for  each  of  the  two  modes,  corresponding 
respectively  to  boxes  S1  and  Sq  in  Fig.  4. 

The  step  junction  is  a  lossless  asymmetric  discontinuity,  and  it 
therefore  requires  three  real  quantities  for  its  characterization.  A 
useful  equivalent  circuit  representation  for  that  type  of  structure  is 
the  one  shown  in  Fig.  5.  It  has  been  found  by  experience  with  careful 
measurements  on  a  variety  of  3tep  junction  discontinuities  in  rectangu¬ 
lar  waveguide  [14]  that  the  series  reactance  X  is  always  very  small, 

9 

and  that  for  most  situations  it  may  be  safely  neglected.  The  represen¬ 
tation  in  Fig.  5(b)  thus  conveniently  reduces  to  a  shunt  network 
comprised  of  a  shunt  susceptance  B  and  a  transformer  with  turns  ratio  n. 


-  12  _ 


The  step  junction  and  a  rigorous  equivalent  circuit  representation  for 
For  most  steps,  Xg  is  very  small  and  can  be  neglected. 


(a)  The  Transformer  Turns  Ratio 

From  the  equivalent  circuit  in  Fig.  5(b),  after  setting  X 

9 

see  that  the  transformer  turns  ratio  is  given  simply  by 


where  the  voltages  V  and  V1  are  given  by 

a*/2 

V(0)  -  /  E»(x,0,z)  •  e*(x,z)dxdz 

-a */2 

a'/2 

V'(0)  -  /  EL(x,0,z)  •  £’*(x,z)dxdz 

-a'/2 


-  0,  we 

(12) 

(13) 


if  we  choose  y-0  to  define  the  plane  of  the  step  Junction  shown  in  Fig. 
5(a),  and  if  the  step  is  uniform  along  z.  For  simplicity,  we  choose  the 


aperture  transverse  electric  field  E,.  to  be  the  same  as  the  mode  field 

— t 

in  the  guide  of  narrower  width,  a'. 

For  the  mode  i-1 ,  we  would  then  write 

-JV  IX  *Jkzz 

E. (x,0,z)  -  E  (x)e  x  -  A  sin  e  x  (1*0 

“w  3  “O  a  — O 

On  use  of  (7)  for  the  mode  function  and  relations  (13)  and  (1*1),  expres¬ 
sion  (12)  for  n  simplifies  to 


a’/2  [2  .  *x  . 
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a  V2  a 
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/  IF sin  r 
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.  nx  . 
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so  that 
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where  the  superscript  S,  which  we  have  just  added,  signifies  "step". 

Proceeding  similarly  for  the  i-0  mode,  one  finds,  using  (9)  for  the 
mode  function, 


S  [-aT 
0  “  J  a 


(b)  The  Shunt  Susceptances 

Since  the  fields  for  the  i-1  and  i-0  transverse  modes  are  dif¬ 
ferent,  the  step  junction  shunt  susceptances  will  also  be  different. 
The  normalized  susceptance  for  the  i-1  mode  has  been  derived  pre¬ 
viously  in  reference  15  in  connection  with  the  transverse  equivalent 
network  for  the  dominant  mode  alone.  The  expression  is 


0.55kyi|icot^ 
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The  factor  0.55  arises  because  (17)  is  obtained  from  an  expression  for 
the  normalized  susceptance  of  a  zero- thickness  transverse  symmetrical 
window  of  relative  aperture  a'/a  subjected  to  i-1  mode  incidence.  From 
stored  power  considerations,  presented  in  reference  15,  and  from  prior 
experience  with  careful  measurements  on  step  Junctions  of  various  types 
in  rectangular  waveguide  Cl1*],  it  was  found  that  the  factor  0.55 
represents  a  good  rule  of  thumb  when  extrapolating  from  a  zero-thickness 
structure  to  a  step  junction  of  similar  dimensions.  The  factor  ky1 
enters  because  it  is  the  wavenumber  of  the  transmission  line  in  the  y 
direction  for  the  i=l  mode. 


For  the  i»0  mode  in  the  y  direction,  we  can  make  use  of  expressions 
for  the  normalized  susceptance  available  in  the  Waveguide  Handbook  [16]. 
One  may  use  either  an  expression  in  Sec.  5.26  on  pages  307  and  308, 
which  applies  directly  to  a  step  junction,  or  expression  2(a)  of  Sec. 
5.1  on  page  218,  after  multiplication  by  0.55  for  the  reasons  stated 
above.  For  convenience,  we  employ  a  simplified  version  of  the  latter 
expression.  Our  result  for  the  normalized  susceptance  of  the  step  junc¬ 
tion  when  the  i«0  mode  is  incident  is 


-  0.55  k  —  in  esc  ~ 
Yq  yO  ir  2a 


Expressions  (17)  and  (18)  are  seen  to  resemble  each  other  somewhat, 
but  they  have  different  functional  dependences  with  respect  to  a'/a.  In 
the  range  of  small  steps,  however,  they  are  really  more  similar  than 
they  appear.  They  may  be  rewritten  in  the  forms 

0.55  ky1  f  ta„2[f  O  -f>]  09) 

^■°-55k,orl"«0[f  <’  -f>]  (2°> 

which  become,  in  the  limit  of  a'/a-O  (small  steps), 
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Both  susceptances  are  seen  to  approach  zero  as  the  square  of  (a-a'),  in 
agreement  with  what  one  obtains  from  small  obstacle  considerations. 


3.  The  Asymmetrical  Coupling  Strip 

We  next  derive  the  equivalent  network  for  the  asymmetrical  strip  of 
width  6  (see  Fig.  3)  which  couples  the  exciting  i-1  transverse  mode  to 
the  radiating  i-0  transverse  mode.  We  should  also  recall  that  the  strip 
is  located  in  the  narrow  parallel  plate  region  of  width  a',  where  the 
i=1  transverse  mode  is  below  cutoff,  so  that  the  strip  is  excited  by  an 
evanescent  field.  Thus,  the  strip  couples  a  propagating  mode  (i*0)  and 
an  evanescent  mode  (i«1). 

A  network  which  accurately  describes  such  coupling  for  arbitrary 
values  of  6/a  might  be  rather  involved,  but  we  seek  in  our  leaky  wave 
antenna  to  produce  small  leakage  per  unit  length,  so  that  the  radiated 
beams  can  be  narrow.  As  a  result,  we  are  concerned  primarily  with  small 
values  of  6/a.  It  is  appropriate  therefore  to  employ  a  small  obstacle 
approach,  which  results  in  a  simple  network  representation  for  the  cou¬ 
pling  effects.  The  small  obstacle  formulation  must  be  valid,  however, 
in  a  multimode  context,  where  one  of  the  modes  is  below  cutoff. 

Fortunately,  a  formulation  which  supplies  most  of  the  requirements 
is  available  in  the  literature  [17].  In  addition,  a  very  simple  and 
convenient  network  form  is  presented  that  is  valid  for  transverse  obsta¬ 
cles,  which  is  precisely  our  case.  The  requirement  that  is  lacking  is 
the  condition  that  one  of  the  modes  is  below  cutoff,  but  the  formulation 
can  be  modified  easily  to  account  for  that  added  step. 

The  small  obstacle  approach  considers  that  the  fields  incident  on 
the  small  obstacle  redistribute  the  surface  charges  such  that  equivalent 
electric  and  magnetic  dipoles  are  created  whose  strengths  are  propor¬ 
tional  to  the  electric  and  magnetic  polarizabilities  of  the  obstacles. 
This  approach  assumes  that  these  polarizabilities  can  be  computed  under 
"static"  conditions  with  the  obstacle  in  "free  space,"  meaning  not  near 
to  any  walls.  The  expressions  for  the  polarizabilities  are,  as  a 
result,  very  simple  in  form.  One  then  computes  the  fields  excited 
within  the  waveguide  by  the  equivalent  dipoles  mentioned  above,  so  that 
the  waveguide  mode  functions  enter  into  the  calculations.  The  field 
components  that  are  taken  into  account  in  determining  the  influence  of 
the  obstacle  are  the  tangential  component (s)  of  the  electric  field  and 
the  normal  component  of  the  magnetic  field.  For  the  asymmetrical  strip 


of  width  6,  the  field  components  are  E  and  H  . 

x  y 

Small  obstacle  theory,  as  indicated  above,  assumes  that  the  obsta¬ 
cle  is  far  from  the  walls,  but  in  our  problem  the  coupling  strip  is 
located  at^  a  wall.  It  is  necessary  to  employ  symmetry  to  modify  the 
cross  section  so  as  to  satisfy  this  consideration.  As  seen  in  Fig.  6, 
the  asymmetrical  strip  of  width  6  in  a  waveguide  of  width  a'  with  the 
i=1  mode  incident  is  equivalent  to  a  centered  symmetrical  strip  of  width 
26  in  a  waveguide  of  width  2a'  with  the  i-2  mode  incident.  The  calcula¬ 
tion  is  therefore  made  first  for  the  latter  set  of  specifications,  and 
then  the  equivalence  is  used  to  adapt  the  result  to  the  guide  of  width 


Reference  17  shows  that  the  coupling  network  can  take  on  the  simple 
form  appearing  in  Fig.  7  since  the  small  obstacle  is  transverse.  That 
reference  also  presents,  in  its  Eq.  (37),  a  simple  approximate  expres¬ 
sion  for  the  admittance  element  that  couples  two  propagating  modes 
in  a  multimode  waveguide,  namely, 

yij  ■  \io  \jo  *  E  <23> 

In  this  expression,  we  have  changed  z  to  y,  in  conformity  with  our  coor¬ 
dinate  system  in  Fig.  3.  and  we  have  added  the  factor  w  which  reference 
17  inadvertently  omitted.  The  Y^  and  Yj  represent  the  characteristic 
admittances  for  the  two  modes,  m  and  P  are  the  magnetic  scalar  and 
electric  tensor  polarizabilities,  and  h  and  e  represent  the  magnetic  and 
electric  mode  functions.  The  subscript  o  means  that  the  value  of  the 
mode  function  is  taken  at  the  center  of  the  obstacle,  and  the  asterisk 
signifies  complex  conjugate.  The  particularly  simple  form  in  Fig.  7 
means  that  y^  can  be  written  as 

yU  ■  yoo  "l  "j  ^ 

where  the  superscript  0  designates  a  turns  ratio  representative  of  the 
"obstacle,"  as  distinguished  from  those  for  the  step  junctions.  (It 
should  be  mentioned  that  in  Fig.  5  of  reference  17  the  turns  ratios 
should  be  specified  as  1 :n  rather  than  n:1.) 


When  one  of  the  modes  that  are  coupled  together  is  below  cutoff, 
expression  (23)  must  be  modified  by  writing  Y^  rather  than  Y^  if  the 
ith  mode  is  the  one  below  cutoff. 


Fig.  6  Left:  Actual  guide  with  asymmetrical  strip  of  width  6.  Right:  Equivalent 
double  width  guide  with  centered  symmetrical  strip  of  width  26. 


Fig.  7.  Simple  equivalent  network,  based  on  a  small  obstacle  approach,  for  the 
asymmetrical  strip  that  couples  the  bound  i  =  1  and  the  radiating  i  =  0 
transverse  modes. 


For  the  asymmetrical  strip,  the  field  components  required  are  Ex  and  H  , 
so  that  expression  (23)  reduces  to 

hj  -  Vj\1o»>Jo  *  J«CP  exioexjo  <25) 

where  -  p  x ,  and  we  have  added  the  complex  conjugate  designation 
on  Yj.  The  polarizabilities  for  a  long  strip  of  width  d  are 

p  -  -m  -  it(d/2)2  (26) 

per  unit  length  of  strip,  in  the  notation  of  reference  17,  and  used 
above.  From  Fig.  6(b)  we  see  that  d/2  -  6,  so  that 

p  -  -ra  -  (27) 


The  mode  functions  for  the  i-2  mode  in  the  waveguide  of  width  2a', 
as  in  Fig.  6,  are  obtained  in  a  manner  identical  to  those  given  above 
for  i-1 ,  in  Eqs.  (7)  and  (8),  subject  to  normalization  (5).  They  are 
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which  are  then  taken  at  x-0  when  used  in  (25).  Similarly,  the  mode 
functions  for  the  i*0  mode  in  the  waveguide  of  width  2a  become 
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The  values  of  kz  and  ky^  above  are  the  same  as  those  for  the  waveguide 

of  width  a',  while  ky2  for  the  guide  of  width  2a'  is  the  same  as  ky1  for 

the  guide  of  width  a'.  In  this  problem,  k  and  k  -  are  real  and  k  ,  is 

z  yu  yd 

imaginary. 

Since  ky2  and  Y2  in  the  guide  of  width  2a'  are  equal  to  ky1  and  Y1 
in  the  original  guide  of  width  a',  and  since  ky(j  and  YQ  are  the  same  in 
both  guides,  we  may  Insert  the  coupling  network  in  Fig.  7  directly  into 
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the  transverse  equivalent  network  in  Fig.  4  after  identifying  the 
appropriate  equivalences  in  the  results  obtained  from  Eq.  (25).  In 
using  (25),  however,  we  shall  assume  that  the  two  modes  involved  are  the 
i«0  and  i-2  modes  in  the  guide  of  width  2a';  subsequently,  the 
equivalences  will  be  made  to  adapt  the  results  to  the  original  i-0  and 
i « 1  modes . 

Although  either  y  or  y22  can  be  u3ec*  *n  the  coupling  network,  we 
have  chosen  yQQ.  When  the  mode  functions  in  (30)  and  (31),  and  relations 
(11)  and  (27),  are  inserted  into  (25),  one  obtains 
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The  turns  ratio  n^  is  obtainable  from  (2*0  as 
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so  that  y22  must  also  be  derived.  Using  (28),  (29),  (2)  and  (27)  in 

(25),  we  find 
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which  differs  from  yQ0  in  (32)  by  only  a  factor  of  2.  From  (33).  (32) 

and  (3*0,  we  thus  deduce  that 


nr 


In  view  of  the  equivalences  mentioned  earlier,  and  the  forms  of  yQQ 
and  n2,  we  can  now  identify  y^^  and  n2  as  computed  above  with  the  yQQ 
and  n^  required  in  Fig.  7,  and  directly  employ  relations  (32)  and  (35). 

The  simple  form  of  Fig.  7  is  a  familiar  one  in  equivalent  network 
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formulations,  and  it  requires  that  the  condition  y00  y22  ~  Yq2  be 
satisfied.  On  use  of  (25),  one  finds  readily  that  the  condition  indeed 
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applies  here, 
is  given  by 


That  formulation  also  indicates  that  the  turns  ratio  n„ 
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Again,  from  (25)  and  (36)  one  readily  verifies  that  (35)  is  correct. 
The  coupling  network  in  Fig.  7  is  therefore  self-consistent  and  valid  in 
the  small-obstacle  limit. 


4 .  The  Radiating  Open  End 

In  the  narrower  regions  of  the  groove  guide,  of  width  a',  the  i  =  1 
transverse  mode  is  evanescent  and  the  i-0  transverse  mode  is  propagat¬ 
ing.  The  length  e",  shown  in  Fig.  3,  is  assumed  to  be  sufficiently  long 
that  the  fields  in  the  evanescent  i-1  mode  have  decayed  to  negligible 
values  as  they  reach  the  open  end.  In  the  transverse  equivalent  net¬ 
work,  whose  form  is  presented  in  Fig.  4,  the  i-1  transmission  line  is 
therefore  taken  to  be  infinite,  and  a  box  labeled  YR  is  placed  at  the 
end  of  the  i»0  transmission  line  to  represent  the  admittance  effect  of 
the  open  end  on  that  propagating  mode. 


For  box  YR  in  Fig.  4,  we  require  the  equivalent  circuit  for  a  TE 
mode  (which  is  what  the  H  z^-type  i-0  mode  reduces  to)  incident  in  the  y 
direction  on  an  open-ended  parallel  plate  guide  of  width  a',  with  the 
electric  field  perpendicular  to  the  plates.  Such  an  equivalent  circuit 
is  already  available  in  the  Waveguide  Handbook  [16],  on  page  179,  in  the 
form  of  a  susceptance  and  a  conductance  in  parallel  across  the  end  of 
the  transmission  line.  Our  notation  differs  in  minor  ways  from  that  in 
the  Waveguide.  Handbook,  but  we  employ  Eqs.  (3)  and  (4),  together  with 
(1)  and  part  of  (2),  on  page  179.  The  biggest  difference  is  that  their 
A'  becomes  our  2ir/k^Q.  Our  admittance,  comprised  of  CR  +  JBR,  where  the 
subscript  R  signifies  "radiating,"  is  taken  directly  at  the  open  end, 
corresponding  to  their  reference  plane  T'. 


In  our  notation,  the  relevant  relations  become 
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cosh  (k^0a'/2)  +  cos  2k^Qd 


where 


2k'  d  - 
yO 


k'  ^ 
yO  it 


An 


_  ?S  ( k  *  ——  *0  0) 

a'  i  'yO  2  it  fV,U> 


Yk 1 

T  yo 


(37) 

(38) 


(39) 
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with  e  =  2.718,  Y-1.781,  and 


S^tjO.O)  -  I  (sin-1 
n-1 


(40) 


The  above  relations  have  been  obtained  by  the  transform  method,  and 
are  rigorous  in  the  range  k^0a'/2ir  <1.  We  have  found,  however,  that 
the  values  of  k  -  8-ja  computed  via  the  transverse  equivalent  network 

Z 

are  affected  negligibly  when  the  infinite  sum  in  (40)  is  replaced  by 
its  first  term  only.  We  have  therefore  replaced  (40)  by 

S^tjO.O)  =  sin-1t  -  t  (41) 

in  our  dispersion  relation. 

5.  The  Complete  Transverse  Equivalent  Network 

In  subsections  1  through  4  above,  we  have  presented  expressions  for 
all  the  constituents  of  the  transverse  equivalent  network  whose  basic 
form  appeared  in  Fig.  4.  The  final  form  of  that  network  is  now 
presented  in  Fig.  8. 

The  two  transmission  lines  comprising  the  network  in  Fig.  8 

represent  the  transverse  i»1  mode,  which  is  the  exciting  mode,  and  the 

transverse  i=0  mode,  which  is  the  mode  created  by  the  added  coupling 

strip  and  which  produces  the  radiation.  Both  transverse  modes  propagate 

in  the  y  direction,  and  they  combine  to  produce  the  net  longitudinal  TE 

mode  in  the  z  direction,  with  propagation  constant  k  the  value  of 

z 

a  representing  the  leakage  per  unit  length  from  the  radiating  open  end. 
As  mentioned  earlier,  the  transverse  equivalent  network  in  Fig.  8  is 
descriptive  of  the  cross  section  of  the  antenna  structure  in  Fig.  3,  but 
it  is  placed  on  its  side  for  clarity. 

The  propagation  wavenumbers  in  the  two  transverse  transmission 
lines  are  related  to  k  by 


Wavenumbers  k  ,  k  n  and  k'  are  basically  real,  and  k',  is  basi- 
yl  ,  yO  yO  yl 

cally  imaginary.  The  characteristic  admittances  are  related  to  the 
respective  wavenumbers  by  (2),  although  those  for  the  i-0  line  can 
employ  the  simpler  form  (11)  in  view  of  (44). 

Boxes  S1  and  SQ  in  Fig.  4,  representing  the  effect  of  the 
step  junctions  on  the  i-1  and  i-0  modes,  respectively,  take  the  form  in 
Fig.  8  of  shunt  susceptances  and  BQ  together  with  ideal  transformers 
of  turns  ratios  and  nQ.  Expressions  for  the  shunt  susceptances  B1/Y1 
and  Bq/Yq,  normalized  to  their  respective  characteristic  admittances, 
are  given  by  (17)  and  (18),  or  in  alternative  forms  (19)  and  (20).  The 
turns  ratios  appear  in  (15)  and  (16). 

The  terminal  admittance  CR  ♦  jBR,  representing  the  effect  of  the 
open  end  on  the  radiating  i-0  transmission  line,  is  described  by  rela¬ 
tions  (37)  and  (38),  together  with  (39)  and  (41). 

The  influence  of  the  added  coupling  strip  in  Fig.  3  is  expressed  in 

basic  form  by  the  box  labeled  Ob,  representing  "obstacle,"  in  Fig.  4, 

and  more  explicitly  by  the  simple  network  in  Fig.  7.  The  network  in 

Fig.  7  is  then  employed  as  part  of  the  overall  network  in  Fig.  8. 

Expressions  for  the  parameters  yrtr.  and  n?  of  this  coupling  network  are 

uu  i  0  0 

presented  in  (32)  and  (35);  the  discussion  there  explains  why  n1  and  n2 

are  the  same  in  the  context  of  (35). 

We  are  pleased  that  we  were  able  to  obtain  all  of  the  constituents 
of  the  transverse  equivalent  network  in  Fig.  8  in  closed  form.  Despite 
the  simple  closed-form  nature  of  the  expressions  for  the  step  junction 
and  the  radiating  open  end,  those  expressions  are  believed  to  be  suffi¬ 
ciently  accurate  to  yield  reliable  numerical  results  for  the  propagation 
behavior  of  the  leaky  groove  guide  antenna.  The  coupling  network 
representing  the  added  asymmetrical  strip  was  derived  using  small  obsta¬ 
cle  theory,  and  was  found  to  be  self-consistent  and  valid  for  small 
values  of  strip  width  6.  The  results  for  0  -  ja  for  small  strip  widths 
should  therefore  be  quite  accurate,  but  the  leakage  rate  o  would  be 
small  for  narrow  strips,  resulting  in  a  long  antenna  and  therefore  a 
rather  narrow  radiated  beam.  From  other  considerations,  which  are 
touched  upon  in  Section  E,  we  believe  that  for  wider  strips  the  leakage 
rate  is  actually  greater  than  that  computed  using  the  small  obstacle 
approach.  In  other  words,  the  actual  strip  width  required  for  a  wider 
beam  is  somewhat  smaller  than  that  indicated  by  the  theory,  meaning  that 
the  range  of  practical  beam  widths  is  greater  than  that  predicted  by  the 


D.  THE  DISPERSION  RELATION 


1 .  The  Dispersion  Equation 

The  dispersion  relation  yields,  as  a  function  of  frequency  and  the 
dimensional  parameters,  the  behavior  of  the  phase  constant  B  and  the 
attenuation  constant  a  of  the  longitudinal  TE  mode  that  propagates  along 
the  z  direction  in  Fig.  3-  From  this  information  the  leaky  wave  antenna 
can  be  designed  according  to  prescribed  specifications.  The  relations 
among  the  frequency  (or  free  space  wavenumber  k,  since  k  -  2irf/c),  the 
longitudinal  wavenumber  k^  *  B_ja,  and  the  transverse  wavenumbers  are 
given  in  ( ^2 )  to  (44).  The  transverse  wavenumbers  ky1 ,  k'y1  and 
kyQ  (=k'yQ)  are  related  by  a  dispersion  equation  that  is  readily  deriv¬ 
able  from  the  transverse  equivalent  network  in  Fig.  8  upon  application 
of  the  transverse  resonance  condition 

7(T)  +  7(T)  -  0  (45) 

Since  relation  (45)  represents  a  free  resonance,  reference  T  may  be 
chosen  anywhere  in  the  network;  our  choice  is  indicated  in  Fig.  8. 

The  application  of  (45)  to  the  transverse  equivalent  network  in 
Fig.  8  yields  the  dispersion  equation  (46)  shown  on  the  next  page. 

The  dispersion  behavior  is  then  obtained  from  this  equation  and 
those  in  (42)  to  (44),  also  recognizing  that  Y'0  -  Yq  . 

2.  Numerical  Procedure 


By  following  the  numerical  procedure  given  by  Bardati  and  Lam- 
par  iello  [18],  the  roots  of  the  complete  transcendental  dispersion  equa¬ 
tion  (46)  on  the  complex  plane  have  been  found  when  parameters  such  as 
the  frequency  f,  the  width  b  of  the  groove,  the  length  c"  of  the  upper 
part  of  the  waveguide,  and  others,  are  varied. 

The  whole  equation  (46)  can  be  regarded  as  an  implicit  function  of 
two  variables  set  equal  to  zero: 

W(z,p)  -  0 

where  z  is  the  actual  unknown  k^b,  and  P  stands  for  the  parameter  which 
is  to  be  varied.  The  problem  of  finding  the  roots  of  this  equation  is 


that  of  making  explicit  the  implicit  function  as  a  function  z  of  the 
parameter  p: 

z  -  z(p) 

If  one  root  of  the  equation,  say  zQ  for  a  given  value  of  the  parameter 

pQ,  is  known  on  the  basis  of  the  implicit  function  theorem  [19],  another 

root  z  exists  inside  a  circle  with  center  at  z  for  a  different  value  of 

o 

the  parameter  p  -  pQ  +  Ap.  The  circle  should  be  sufficiently  small  so 
as  to  contain  only  one  root  for  each  of  the  values  pQ  and  pQ  +  Ap  of  the 
parameter. 

Using  this  procedure,  let  us  assume  that  we  know  the  zeros  of  the 
transverse  resonance  equation  for  the  i-1  mode  when  the  coupling  strip 
is  absent  and  the  length  c"  of  the  upper  part  of  the  waveguide  is  infin¬ 
ite.  Now,  let  us  write 

W(z,h)  -  A(z)  +  -  0  (H7) 

where  A(z)  is  the  function  that  refers  to  the  above-mentioned  i«1  mode. 
The  zeros  of  this  equation  are  identical  to  those  of  the  complete 
transverse  resonance  equation  when  the  parameter  h  is  set  equal  to 
unity.  Starting  from  the  known  root  of  the  equation  for  h-0,  it  is  pos¬ 
sible  to  calculate  the  final  solution  by  varying  the  parameter  h  from  0 
to  1 .  If  we  consider  the  function  [W(z,h)]h-1  as  a  function  of  any 
other  parameter  which  is  of  interest  to  be  varied,  the  frequency  for 
example,  the  equation 

W(z,f)  -  0 

can  be  solved  starting  from  the  now  known  root  for  a  given  frequency; 
that  is,  the  explicit  function. 

z  -  z(f) 

can  be  numerically  constructed. 

It  is  worthwhile  to  note  that  the  numerical  method  employed  for 
solving  transcendental  equations  [20]  provides  roots  with  very  high  pre¬ 
cision.  For  the  function  arising  here,  the  real  and  the  imaginary  parts 
of  the  function  calculated  at  the  zero  are  at  least  12  orders  of  magni¬ 
tude  lower  than  the  minimum  value  of  the  modulus  of  the  function  along  a 


E.  NUMERICAL  VALUES  AND  DESIGN  CONSIDERATIONS 


1 .  Some  General  Properties  of  Leaky  Wave  Antennas 

In  order  to  place  the  numerical  values  to  be  presented  below  in 
proper  context,  it  is  important  to  recall  a  few  general  properties  of 
leaky  wave  antennas. 

A  leaky  wave  antenna  is  basically  an  open  waveguide  possessing  a 
mechanism  which  permits  a  slow  leakage  of  power  along  the  length  of  the 
waveguide.  The  leaky  wave  that  exists  at  the  opening  along  the  length 
of  the  waveguide  provides  the  aperture  field  distribution  of  the 
antenna,  and  the  length  of  the  waveguide  over  which  the  leaky  wave  con¬ 
tains  significant  power  defines  the  antenna  aperture  length.  It  is  cus¬ 
tomary  to  retain  a  length  for  which  90$  or  95%  of  the  power  can  leak 
away;  the  remaining  10$  or  5$  is  absorbed  in  a  load  placed  at  the  end  of 
the  waveguide. 

The  rate  of  leakage  per  unit  length  along  the  waveguide  is  given  by 
the  attenuation  constant  a  of  the  complex  propagation  constant  k  (»fj-ja) 
of  the  leaking  waveguide  mode.  If  the  value  of  a  is  small,  the  antenna 
aperture  will  be  long,  and  the  far  field  radiation  pattern  will  possess 
a  narrow  beam.  If  a  is  large,  the  power  is  radiated  more  rapidly  along 
the  antenna  length  and  a  wider  radiated  beam  will  be  obtained. 

The  design  of  a  leaky  wave  antenna  proceeds  by  fir3t  specifying  the 
desired  performance  characteristics  of  the  antenna:  the  angle  of  maximum 
radiation,  the  beam  width,  and  the  side  lobe  properties.  The  angle  0m 
that  the  maximum  of  the  main  beam  makes  with  the  normal  to  the  antenna 


aperture  (broadside  direction)  is  given  approximately  by 


sin  0  -  6/k  -  A/A„  (48) 
m  g 

where  B  is  the  phase  constant  of  the  complex  propagation  constant  k  . 
The  beam  width  A6  is  related  linearly  to  the  reciprocal  of  the  aperture 
length  L,  and  a  rule  of  thumb  for  A0  is 


A0 


A 

l  cos  e 

m 


(49) 


Relation  (48)  above  is  summarized  pictorially  in  Fig.  9. 


For  the  side  lobe  properties,  it  is  well  known  that  the  antenna 
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Fig.  9  Pictorial  representation  for  sin  0mS:  g/k,  where  0m  is  the  angle  that 

the  maximum  of  the  radiated  beam  makes  with  the  normal  to  the  antenna 
aperture  surface,  and  g  is  the  phase  constant  of  the  leaky  wave  that 
travels  along  the  aperture  surface. 


aperture  distribution  must  be  tapered  in  a  specified  fashion  in  order  to 
achieve  the  desired  side  lobe  requirements.  The  necessary  tapers 
corresponding  to  various  side  lobe  specifications  appear  in  standard 
antenna  texts.  When  the  aperture  distribution  is  tapered,  relation  (49) 
for  the  beam  width  is  affected  somewhat  numerically,  but  not  more  than 
±25$  depending  on  the  taper,  so  that  (49)  remains  a  useful  and  reason¬ 
ably  accurate  indication  of  the  beam  width. 

When  the  geometry  of  the  waveguide  is  uniform  along  its  length,  the 

aperture  distribution  consists  of  a  slow  exponential  decay.  To  produce 

the  taper  needed  for  the  side  lobe  requirements,  it  is  necessary  to  vary 

the  a  of  the  waveguide  along  the  length  L,  and  therefore  to  vary  the 

transverse  dimensions  in  some  fashion  along  the  waveguide  length.  The 

phase  constant  8  is  maintained  constant  along  the  length  so  that  all 

parts  of  the  antenna  aperture  radiate  at  the  same  angle  6  .  The  rela- 

m 

tion  between  a(z),  where  z  is  the  longitudinal  coordinate,  and  the 
desired  taper  along  z  of  the  antenna  aperture  distribution  is  also 


available  in  the  literature. 


The  dispersion  relation  (46)  for  the  leaky  groove  guide  yields  the 
values  of  a  and  8  corresponding  to  specific  geometrical  parameters  in 
the  cross  section.  That  relation  must  be  used  in  a  specific  design  to 
produce  the  a(z)  needed  while  maintaining  B(z)  constant.  Part  of  the 
design  considerations  to  be  presented  below  are  aimed  at  recommending 
which  geometrical  parameter  or  parameters  should  be  varied  to  change  a 
but  not  8. 


The  remainder  of  the  report  discusses  the  behavior  of  the  leaky 
groove  guide  when  the  cross  section  is  maintained  constant,  and  it 
stresses  the  variation  of  8  and  a  with  frequency  and  with  various  dimen¬ 
sional  parameters.  If  we  assume,  as  is  common,  that  the  antenna  aper¬ 
ture  of  length  L  radiates  90 ?  of  the  power  (with  the  remaining  10? 
absorbed  by  a  load),  then  the  power  level  at  the  end  of  the  antenna 
aperture  is  10  dB  down  from  its  value  at  the  beginning.  In  the  numeri¬ 
cal  values  that  follow,  we  shall  present  the  a  value  in  the  form  ot/k  in 
nepers,  so  that  the  value  of  a  in  dB  per  wavelength  is  2u(8.686) (a/k)  ** 
54.6  (a/k).  The  antenna  length  L  in  free  space  wavelengths  is  then 

approximately 


L  _  10  1  -  0,18 

A  "  54.6  a/k  a/k 


(50) 


The  precise  value  of  L  will  depend  on  how  the  aperture  distribution  is 
tapered.  If  the  percentage  of  power  radiated  is  different  from  90?,  the 
length  L  in  (50)  will  differ  somewhat  numerically,  but  may  be  readily 
computed. 


In  summary,  the  values  of  8/k  and  a/k  that  are  presented  below  per¬ 
mit  us  to  determine  readily  the  values  of  0_  from  (48),  the  aperture 

•  m 

length  L  from  (50),  and  the  beamwidth  66  from  (49). 


2.  Optimization  and  Other  Design  Considerations 

In  order  to  systematically  design  radiation  patterns,  one  must  be 
able,  as  explained  in  the  section  Just  above,  to  taper  the  antenna  aper¬ 
ture  amplitude  distribution  while  maintaining  the  phase  linear  along  the 
aperture  length,  i.e.,  one  must  be  able  to  vary  a  while  keeping  e  the 
same.  Fortunately,  several  parameters  can  be  varied  that  will  change  a 
while  affecting  6  hardly  at  all;  those  parameters  are  6,  c,  c'  and  c", 
as  is  evident  from  the  numerical  values  given  in  the  next  section. 


Practical  considerations  may  rule  out  parameters  6  and  c’;  lengths  c  and 
c"  must  be  large  enough  that  changing  them  does  not  affect  the  evanes¬ 
cent  field  distribution  of  the  i«l  transverse  mode. 

In  the  design,  one  first  chooses  the  width  a  and  adjusts  a'  and  b 
to  achieve  the  desired  value  of  8/k,  which  is  determined  essentially  by 
the  i»1  transverse  mode.  That  value  of  8/k  immediately  specifies  the 
angle  of  the  radiated  beam.  It  is  then  recognized  that  the  value  of  a 
can  be  increased  if  the  coupling  strip  width  6  is  increased,  or  if  the 
distance  c'  between  the  step  junction  and  the  coupling  strip  is 
decreased,  since  the  coupling  strip  is  excited  by  the  i=1  transverse 
mode,  which  is  evanescent  away  from  the  step  junction  in  the  outer 
regions. 

The  selection  of  the  values  of  the  remaining  dimensional  parame¬ 
ters,  c  and  c" ,  requires  special  considerations.  Since  the  i*0 
transverse  mode  is  above  cutoff  in  both  the  central  and  outer  regions  of 
the  guide  cross-section,  a  standing  wave  effect  is  present  in  the  i=0 
transmission  line.  As  a  result,  a  short  circuit  can  occur  in  that 
transmission  line  at  the  position  of  the  coupling  strip  of  width  6,  and 
the  value  of  a  then  becomes  zero.  Hence,  we  must  choose  the  dimensions 
to  avoid  that  condition,  and  in  fact  to  optimize  the  value  of  a.  There¬ 
fore,  after  c'  is  specified,  the  length  c  must  be  determined  such  that 
the  standing-wave  effect  mentioned  above  optimizes  the  value  of  a.  If  c 
is  sufficiently  long,  it  will  affect  only  the  i-0  transmission  line  and 
influence  8  negligibly.  The  length  cw  also  affects  a  strongly  and  6 
weakly,  and  it  also  must  be  optimized  because  another,  although  milder, 
standing  wave  exists  between  the  coupling  strip  and  the  radiating  open 
end. 


It  is  important  to  realize  that  the  dimensions  for  optimization  are 
Independent  of  frequency,  since  the  transverse  wavenumbers  are  all  fre¬ 
quency  independent.  Of  course,  when  the  frequency  is  altered  the  values 
of  8  and  o  will  change,  but  the  dimensional  optimization  remains  undis¬ 
turbed. 

The  frequency-independent  nature  of  the  dimensional  optimization 
can  be  put  to  practical  advantage.  Suppose  that  a  curve  of  a/k  vs.  c' 
and  the  corresponding  curve  of  a/k  vs.  c,  where  the  values  of  c  are 
optimized,  have  been  obtained  for  a  given  frequency  f ^ .  If  now  similar 
curves  are  desired  for  a  second  frequency  f2»  it  is  not  necessary  to 
repeat  the  optimization  calculations.  One  simply  takes  a  curve  of  a/k 
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vs.  frequency,  and  finds  the  ratio  R  of  a/k  at  each  of  the  two  frequen¬ 
cies  f2  and  f 1 .  The  new  values  of  a/k  at  frequency  f 2,  as  a  function  of 
c'  and  c,  are  found  simply  by  taking  the  values  of  a/k  for  f1  at  each  of 
the  c'  and  c  values  and  multiplying  each  of  them  by  the  ratio  R. 

The  frequency-independent  property  of  the  transverse  wavenumbers 
leads  to  another  simple  relation  which  may  sometimes  be  of  value.  By 
taking  the  imaginary  parts  of  the  wavenumber  relation 

ag  -  Im(ky1 )Re(ky1 )  (51) 


where  k  =  g-ja  and  k  . 

z  yi 


Re(kyl )  ♦  j  Im(ky1 ),  or 

a  C _ 1_ 

k  “  ,2  g/k 


where  C  is  a  factor  independent  of  frequency  since  ky1  is  frequency 
independent.  Relation  (52)  may  be  useful  if  the  values  of  a/k  and  g/k 
are  known  at  a  given  frequency,  and  the  value  of  a/k  is  desired  at 
another  frequency.  Accurate  values  of  g/k  are  obtainable  from  the  com¬ 
plete  dispersion  relation  (46),  but  approximate  values  can  be  found  from 
the  much  simpler  expression 


cot  k 

yi  2 


i  kyi  + 

tnf)2  l1 


where  n1  and  B.|/Y1  are  given  by  (15)  and  (17);  (53)  is  obtained  by 
removing  the  coupling  strip  and  making  c  and  c’  infinitely  long,  so  that 
kyl  and  |k^  |  become  purely  real  here.  An  approximate  value  of  a  for  a 
new  frequency  f2  can  then  be  determined  using  relation  (52),  where  the 
value  of  C  is  obtained  from  the  a/k  and  g/k  values  at  the  given  fre¬ 
quency  f 1 .  Alternatively,  (52)  may  be  rephrased  as 


(a/k). 


2  f~ 


(6/k), 


(a/k), 


A  third  consequence  of  the  frequency  independence  of  the  transverse 
wavenumbers  is  that  the  beam  width  of  the  radiation  remains  essentially 
constant  as  the  beam  is  scanned  by  varying  the  frequency.  Since 


-  3H  - 


using  ( 48) ,  and 


cos20  -  1-sin2e  *  1-(B/k)2 
m  in 


1  —  ( 8/k ) 2  -  (ky1/k)2  +  (n/ka)2  -  (k^/k)2 


using  (42),  we  may  rewrite  (49)  as 


A6-AJL.  (55) 

L  kt1 

where  k  is  the  transverse  wavenumber  corresponding  to  the  i=1 

I 

transverse  mode.  In  the  rephrasing  of  (42)  above,  the  value  of  o  is  set 
equal  to  zero,  so  that  ky1  ,  and  therefore  kfc1 ,  become  purely  real. 

Since  k  is  frequency  independent,  and  k  •=  2tt/A,  we  see  that  the  beam 

Lr  I 

width  becomes 


which  is  independent  of  frequency. 

3.  Numerical  Results 

Using  the  numerical  procedure  described  in  Section  D,2  to  obtain 
numerical  values  from  the  dispersion  relation  (46),  we  have  obtained 
numerical  values  in  graphical  form  of  the  variation  of  a/k  and  B/k  with 
each  of  the  dimensional  parameters.  We  report  here  a  set  of  curves  for 
a  typical  set  of  dimensions. 

Consistent  with  the  discussion  in  the  preceding  section  on  design 
considerations,  we  begin  with  the  width  ratio  a'/a  -  0.68  and  a  height 
b/a  -  0.40  for  the  central  region  of  the  groove  guide  (see  Fig.  3). 
These  two  dimensional  parameters  are  the  principal  ones  that  determine 
the  value  of  B/k.  To  obtain  a  reasonably  large  value  of  a/k,  we  choose 
a  strip  width  6/a  -  0.25  and  its  location  c'/a  -  0.20.  Then,  for  a 
relative  wavelength  A/a  -  1.20,  we  found  that  the  length  c/a  had  to  be 
1.28  for  optimization.  The  corresponding  optimized  value  of  c"/a  turned 
out  to  be  1.50.  These  optimization  considerations  are  discussed  In  sub¬ 
section  2. 

Curves  of  B/k  and  a/k  as  a  function  of  fa  (where  f  is  the  frequency 
in  GHz  and  a  is  the  width  of  the  central  region),  corresponding  to  the 
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above  set.  of  normalized  dimensions,  are  presented  in  Figs.  10(a)  and 
(b).  Figure  10(a)  is  replotted  in  Fig.  11  with  £/k  replaced  by  Bm,  the 
angle  of  the  maximum  of  the  radiated  beam  (see  Fig.  9).  The  relative 
wavelength  A/a  «  1.20,  mentioned  above,  corresponds  to  fa  -  25  in  these 
figures.  It  is  seen  that  the  £  and  a  values  vary  over  a  wide  range, 
with  particular  sensitivity  near  cutoff.  For  these  dimensions,  the 
radiated  beam  can  be  scanned  with  frequency  from  about  15°  to  nearly  60° 
from  the  normal  before  the  next  mode  begins  to  propagate. 

In  order  to  determine  the  variation  of  $/k  and  a/k  with  the  dimen¬ 
sional  parameters,  we  must  select  a  specific  value  of  fa,  or  A/a.  The 
value  A/a  »  1.20  mentioned  earlier  (corresponding  to  fa  -  25)  is  seen 
from  Fig.  10(b)  to  yield  a  reasonably  large  value  of  a/k,  and  yet  not  be 
too  close  to  cutoff.  For  guide  width  a  -  5.0  mm,  this  relative 
wavelength  corresponds  to  a  frequency  of  50  GHz. 

In  Fig.  12(a),  we  present  a  curve  of  a/k  vs.  c'/a,  the  normalized 
distance  of  the  asymmetric  strip  from  the  step  junction,  corresponding 
to  the  relative  wavelength  A/a  -  1.20.  As  expected,  small  changes  in  c' 
produce  large  variations  in  a.  It  is  next  necessary  to  find  the  values 
of  c/a,  the  normalized  length  of  guide  between  the  lower  step  junction 
and  the  short-circuited  bottom  of  the  groove  guide,  in  accordance  with 
the  optimization  discussed  under  2.  The  optimization  curve  of  a/k  vs. 
c/a  is  shown  in  Fig.  12(b).  In  effect,  this  curve  indicates  the  value 
of  c/a  required  once  a/k,  via  c’/a,  is  specified.  It  is  interesting  to 
note  that  the  value  of  c/a  +  c'/a  required  for  optimization  is  nearly 
constant,  so  that  a  curve  of  a/k  vs.  (c+c')/a,  as  in  Fig.  12(c),  may  be 
used  as  an  alternative  to  that  in  Fig.  12(b). 

In  order  to  be  sure  that  higher  mode  coupling  between  the  discon¬ 
tinuities  produced  by  the  asymmetric  coupling  strip  and  the  step  junc¬ 
tion  is  acceptably  low,  we  determined  the  decay  rate  of  the  i«2 
transverse  mode.  Under  the  conditions  discussed  here,  it  was  found  that 
| ky2 I /a  c  0.82,  so  that  the  field  of  that  mode  has  decayed  to  less  than 
0.2  of  its  original  value  over  the  distance  c’/a  -  0.20.  All  the  other 
higher  modes  decay  much  more  rapidly.  The  i-2  mode  contains  only  a 
small  portion  of  the  stored  energy,  however,  so  that  the  coupling  strip 
discontinuity  can  safely  be  regarded  as  isolated  with  respect  to  the 
higher  modes.  The  step  discontinuity,  moreover,  does  not  even  excite 
this  transverse  mode;  the  lowest  higher  mode  permitted  by  its  symmetry 
is  the  i»3  mode,  which  decays  much  more  rapidly.  We  are  thus  safe  in 
using  the  network  in  Fig.  8. 


1.5x10 


10  Variation  of  normalized  phase  constant  p/k  and  leakage  constant  a/k  with  normalized  frequency 
in  the  form  fa,  where  f  is  the  frequency  in  GHz  and  a  is  the  width  of  the  guide's  central  region 
(see  Fig.  3).  The  other  geometrical  parameter  values  are  given  in  the  box  in  the  figure. 


I 


Fig.  1 


i 


aVa=  0.68,  b/a=0.40,  c/a  =1.28 
c'/a  =  0.20,  c'/a*  1.60,  S/a=0.26 


fa  (f  in  GHz) 
(a  in  cm) 


1  Variation  of  the  angle  @m  of  the  maximum  of  the  radiated  beam  with 
the  normalized  frequency  fa,  defined  as  in  Fig.  10. 
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a'/a=0.68,  b/a=0.40,  c*/a=1.50 
S/a=0.25,  X/a=1.20 


(b)  The  corresponding  values  of  c/a  required  to  optimize  the  values  of  ct  'k,  where  c  is  the  length  of  gt 
between  the  lower  step  junction  and  the  bottom  of  the  groove  guide,  fc)  An  alternative  phrasing  for  the 
dimensions  for  optimization,  in  terms  of  the  sum  (c  +  c '  )/a .  The  values  of  the  normalized  wavelength 
and  the  other  dimensions  are  given  in  the  box. 


The  variations  of  8/k  and  a/k  with  the  normalized  groove  guide 
height  b/a  are  shown  in  Fig.  13(a)  and  (b),  where  X/a  is  taken  as  1.20, 
and  the  other  dimensions  are  the  ones  given  earlier,  including  the 
optimization  indicated  in  Figs.  12.  It  is  seen  that  8/k  varies  consid¬ 
erably  with  b/a,  particularly  for  smaller  b/a  values.  The  leakage  con¬ 
stant  a/k  also  varies  strongly  with  b/a,  showing  a  maximum  value  at 
about  b/a  -  0.34,  rather  that  the  0.40  value  taken  above  (and  below)  in 
connection  with  other  figures.  It  is  therefore  possible  to  produce  an 
additional  optimization,  by  changing  the  value  of  b/a  used  in  the  deter¬ 
mination  of  Figs.  12  from  0.40  to  0.34,  but  the  increase  in  a/k  will  not 
be  large.  It  is  also  seen  that  a/k  becomes  zero  in  the  neighborhood  of 
b/a  =  1.4.  The  reason  for  this  zero  is  that  a  null  in  the  standing  wave 
for  the  i=0  transverse  mode,  which  is  above  cutoff  in  both  the  £  and  aj 
regions,  then  occurs  exactly  at  the  location  of  the  asymmetric  strip, 
thereby  shorting  out  the  coupling  between  the  i-0  and  the  i-1  transverse 
modes. 

The  dependences  of  8/k  and  a/k  on  the  normalized  length  6/a  of  the 
asymmetric  coupling  strip  are  shown  in  Fig.  14(a)  and  (b).  As  expected, 
a/k  depends  very  strongly  on  6/a,  whereas  8/k  changes  hardly  at  all. 

The  theoretical  expression  for  the  coupling  strip  wa3  obtained  by 
employing  small  obstacle  theory,  so  that  the  result  becomes  less  accu¬ 
rate  as  6/a  increases  beyond  its  range  of  validity.  Expression  (32)  may 
be  valid  up  to  6/a  *  0.25  or  so.  In  a  practical  configuration,  6/a  may 
of  course  be  made  larger  than  that  value  if  one  wishes  to  obtain  greater 
leakage  per  unit  length,  even  though  our  theory  becomes  inaccurate  in 
that  range.  By  comparison  with  similar  discontinuities,  however,  for 
which  there  are  available  both  small  obstacle  results  and  ones  for  arbi¬ 
trary  dimensions,  one  can  conclude  that  the  small  obstacle  results,  if 
extrapolated  beyong  their  range  of  validity,  yield  values  of  coupling 
which  are  lower  than  the  correct  values.  This  means  that  the  actual  a/k 
values  for  large  6/a  ratios  are  expected  to  be  greater  than  the  values 
computed  using  the  small  obstacle  theory  in  the  range  beyond  its  vali¬ 
dity. 


As  indicated  above,  the  value  of  c/a  must  be  selected  in  accordance 
with  the  optimization  specified  in  Figs.  12.  The  ways  in  which  8/k  and 
a/k  vary  with  c/a  appear  in  Fig.  15(a)  and  (b).  It  is  seen  that  8/k 
varies  very  little  with  c/a,  whereas  a/k  is  sensitive  to  it,  with  nulls 
appearing  at  those  values  of  c/a  for  which  the  coupling  strip  is  shorted 
out. 
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.68,  c/a  =  1.28,  c'/a=0.20 
.50,  b/a  =  0.40,  \/a  =  1.20 


Fig.  14  The  variations  of  g/k  and  a/k  with  the  normalized  length  6/a  of  the 

perturbing  asymmetric  strip;  the  values  of  \/a  and  the  other  dimensions 
are  given  in  the  box. 


The  final  variations  that  we  consider  are  those  of  B/k  and  ct/k  with 
c"/a,  the  normalized  length  between  the  coupling  strip  and  the  upper 
open  end.  Those  variations  are  seen  in  Fig.  16(a)  and  (b).  Because  a 
standing  wave  in  the  i-0  transverse  mode  also  exists  in  that  length,  we 
find  that  o/k  varies  with  c"/a,  so  that  it  is  necessary  to  select  the 
best  value  of  c"/a  to  optimize  the  value  of  a/k.  We  chose  the  peak  at 
c"/a  =  1.50  even  though  it  results  in  a  taller  structure  because  we 
found  that  the  field  of  the  evanescent  i-1  transverse  mode  had  not  yet 
decayed  to  negligible  values  at  the  c"/a  value  corresponding  to  the 
first  peak. 

The  actual  proportions  in  the  cross  section  of  the  groove  guide 
antenna  corresponding  to  the  optimizations  discussed  above  are  shown  in 
Fig.  17.  The  antenna  cross  section  ends  up  being  tall  and  thin  to 
insure  that  the  i=1  transverse  mode,  which  is  evanescent  in  the  regions 
of  width  a.',  has  negligible  amplitude  at  the  very  top  and  bottom. 

In  an  actual  antenna  design,  the  value  of  a  must  be  tapered,  as 
indicated  in  Section  E,1,  but  the  value  of  B  must  remain  essentially 
constant  along  the  antenna  length.  From  the  curves  in  Figs.  13  through 
16,  we  see  that  lengths  c  and  c"  are  suitable  parameters  to  vary  for 
that  purpose. 

For  the  specific  case  corresponding  to  Fig.  17,  the  angle  6m  of  the 
maximum  of  the  radiated  beam  (see  Fig.  9  and  Eq.  (48))  is  about  44°,  the 
length  L  of  the  antenna  is  about  25  free  space  wavelengths,  in  accor¬ 
dance  with  the  conditions  of  Eq.  (50),  and  the  beam  width  is  about  3.2°, 
from  Eq.  (49).  It  should  be  recognized  that  the  optimizations  discussed 
above  are  intended  to  increase  the  value  of  o/k,  and  thus  to  Increase 
the  beam  width.  It  is  easy  to  achieve  a  narrower  beam,  by  simply  relax¬ 
ing  the  optimizations. 

We  have  described  here  a  new  type  of  leaky  wave  antenna  suitable 
for  millimeter  waves.  It  is  based  on  the  groove  guide,  which  is  a  low- 
loss  waveguide,  so  that  the  waveguide  attenuation  will  interfere  negli¬ 
gibly  with  the  leakage  process.  The  structure  is  simple  and  longitudi¬ 
nally  continuous,  rendering  it  reasonable  for  fabrication  at  the  small 
wavelengths  in  the  millimeter-wave  range. 

We  have  also  derived  all  the  constituents  of  an  accurate  transverse 
equivalent  network,  which  yields  a  transverse  resonance  relation  in 
closed  form,  from  which  numerical  design  values  can  be  readily  obtained. 


1.68,  c/a  *  1.28,  c'/a=0.20 
.40  8/a=  0.25,  X/a  =  1.20 


variations  of  f}A  and  a  A  with  the  normalized  length  c#/a  between 
asymmetric  strip  and  the  upper  open  end;  the  values  of  x/a  and  the 
r  dimensions  are  given  in  the  box. 


The  antenna  structure  is  therefore  capable  of  straightforward  under¬ 
standing  and  systematic  design.  It  is  also  sufficiently  flexible  with 
respect  to  the  dimensional  parameters  which  can  be  varied  that  a  reason¬ 
ably  wide  range  of  pointing  angles  and  beam  widths  can  be  achieved  with 
it. 
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